Stabilizing continuous-wave output in semiconductor lasers by time-delayed feedback 



Thomas Dahms, Philipp Hovel, and Eckehard ScholQ 
Institut filr Theoretische Physik, Technische Universitdt Berlin, 10623 Berlin, Germany 

(Dated: October 29, 2008) 

The stabihzation of steady states is studied in a modified Lang-Kobayashi model of a semicon- 
ductor laser. We show that multiple time-delayed feedback, realized by a Fabry-Perot resonator 
coupled to the laser, provides a valuable tool for the suppression of unwanted intensity pulsations, 
and leads to stable continuous-wave operation. The domains of control are calulated in dependence 
on the feedback strength, delay time (cavity round trip time), memory parameter (mirror reflectiv- 
ity), latency time, feedback phase, and bandpass filtering, Due to the optical feedback, multistable 
behavior can also occur in the form of delay-induced intensity pulsations or other modes for certain 
choices of the control parameters. Control may then still be achieved by slowly ramping the injection 
current during turn-on. 



I. INTRODUCTION 

Time-delayed feedback methods have been widely used 
to control unstable dynamics in a variety of different 
fields [1 . In its original form, time-delay auto synchroni- 
sation (TDAS) was introduced by Pyragas [2 to stabilize 
unstable periodic orbits embedded in a chaotic attractor. 
In TDAS, the control force is generated from the differ- 
ence of an output signal s(t) and its counterpart s{t — r) 
some time units r ago. One groundbreaking advantage of 
TDAS is the noninvasiveness of the control, i.e., the con- 
trol force vanishes on target orbit. This control scheme 
was extended by Socolar et a/., by considering all integer 
multiples m of the delay time s{t — mr) — s{t — {m-\-l)r) 
weighted with a memory parameter {extended time- 
delay auto synchronisation^ ETDAS) [3 . This control 
scheme was invented in the context of optical systems like 
semiconductor lasers, where feedback can be realized all- 
optically, for instance, by a Fabry-Perot (FP) resonator, 
which naturally generates an ETDAS control force by 
multiple reflections. Such all-optical noninvasive control 
has indeed been realized experimentally by coupling a 
multisection semiconductor laser with an external Fabry- 
Perot cavity to stabilize unstable steady states [H |5] . 

A simple model describing a semiconductor laser with 
optical feedback from a single mirror was introduced by 
Lang and Kobayashi [6^ . The effect of time-delayed feed- 
back (TDAS) on semiconductor lasers was investigated 
within Lang-Kobayashi (LK) type models [TuSilQi irolITT] . 
as well as within more elaborate device models ^ ^3 
\T3\ [T4[ \TE\ \T6\ [TT] for various configurations, includ- 
ing Michelson interferometers providing a realization of 
TDAS. These findings were supported by experimental 
work [4] |14l [18] . Not only noninvasive control but also 
delay-induced instabilities [I9l [20] and high-dimensional 
chaos resulting from time-delayed feedback were studied 
|2Tj . In particular, feedback-induced stationary external 
cavity modes and their bifurcations in a LK model for 
a laser subject to resonant feedback from a Fabry-Perot 



resonator were treated within the TDAS approximation 

In this work, we will consider a modified LK model 
of a semiconductor laser coupled to a Fabry-Perot res- 
onator (see Fig. [l] for a schematic diagram of the sys- 
tem), and investigate how ETDAS can be successfully 
used to stabilize unstable steady states, i.e., continuous 
wave (cw) emission, of the uncontrolled (uncoupled) sys- 
tem. It thus provides a systematic theoretical framework 
for the type of experimental configuration used in Ref. [4 . 
We calculate the domains of control in dependence on 
various control parameters, latency, and bandpass filter- 
ing, which extends our previous findings from a simple 
linear generic model on time-delayed feedback control of 
a fixed point |22l [23j [24] . In addition, we present numer- 
ical simulations of the full nonlinear model and develop 
strategies for global control of the fixed point even in case 
of delay-induced multistability. 
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FIG. 1: Schematic diagram of a semiconductor laser with 
resonant feedback from a Fabry-Perot resonator. K denotes 
an attenuator, R is related to the mirror reflectivity of the 
external resonator, r is the round trip time of the resonator, 
cp and S are phase shift and latency time due to the distance 
between laser and resonator. 

This paper is organized as follows: In Section [TTj we 
introduce a model of a semiconductor laser with inter- 
nal passive dispersive feedback. This model allows for 
undamped relaxation oscillations for certain choices of 
the parameters. Furthermore, optical feedback from a 
Fabry-Perot resonator is included. In order to inves- 
tigate the properties in the vicinity of the lasing fixed 
point, a linear stability analysis of the system is per- 
formed in Section |III[ both with and without external 
optical feedback. The numerical results are presented in 
Section [TVl and the dependence upon feedback strength. 
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time delay, memory parameter, optical phase shift, and 
latency time is studied. In Section [V| a bandpass filter is 
applied to the feedback signal. Again, numerical results 
show the impact on stability of the lasing fixed point. 
The full nonlinear dynamic behavior of the system is in- 
vestigated in Section [Vl| To overcome the problems of 
delay-induced multistability, simulations of a specifically 
designed turn-on process allowing for successful control 
of the steady state are discussed in Section |VII[ Finally, 
we finish with a conclusion in Section [Villi 



II. THE MODEL 

Semiconductor lasers with external optical feedback 
from a mirror can be described by the Lang-Kobayashi 
model [6]. In dimensionless form, it consists of two dif- 
ferential equations for the slowly varying amplitude (en- 
velope) E{t) of the complex electric field and the reduced 
carrier density (inversion) n{t). 

Here we consider a modification of the LK equations 
appropriate for multisection semiconductor lasers with 
an internal passive dispersive reflector [25 . This is mod- 
eled by a gain function k{n) depending upon the inter- 
nal dispersive feedback from the Bragg grating. Such a 
laser structure allows for more complex dynamic behavior 
including self-sustained relaxation oscillations (intensity 
pulsations) generated by Hopf bifurcations, as has been 
shown in the framework of traveling wave laser models 
[4^ 15 . We are interested in the regime above a supercrit- 
ical Hopf bifurcation where the fixed point in the uncon- 
trolled system is unstable. Combining the rate equation 
for the carrier density from Ref. [25 with the rate equa- 
tion for the complex electric field, we obtain the following 
form of modified LK equations: 

^ = ^il+ia)nE-Eb{t), (la) 

^ = I-n-{l+n)k{n)\E\\ (lb) 

where a denotes the linewidth enhancement factor, / is 
the reduced excess injection current, T is the time scale 
ratio of the carrier lifetime Tc and the photon lifetime 
Tp, and Ei){t) denotes the feedback term, which will be 
described in detail later. 

Here we have scaled (i) time t by the carrier lifetime 
Tc, (ii) carrier density n (in excess of the threshold carrier 
density) by the inverse of the differential gain Gn times 
Tp, and (iii) electric field E by (tcGat)"^^^, so that all 
variables are dimensionless. Moreover, the variables / 
and n are linearly transformed such that both are zero at 
the laser threshold. Note that here, as in Ref. [25], time 
is scaled by Tc rather than by Tp, as often used elsewhere 
(see, e.g. [26]). The function /c(n), which models the 
internal dispersive feedback, is chosen as a Lorentzian, 
as proposed in Ref. [25] : 



where A denotes the height, W is the width, and no 
is the position of the resonance. The parameter ko is 
chosen such that k{0) = 1 at the laser threshold. A 
typical curve k{n) is shown in Fig. [2ja). Througout the 
following we will use the parameters A = 1, W = 0.02, 
and no = -0.034. 
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FIG. 2: (a) Typical shape of the gain function k{n) used in the 
modified LK model Eq. ([T]). Parameters: A = 1, W = 0.02, 
and no = —0.034. (b) Dependence of the damping rate F on 
the choice of the detuning parameter no in the function k(n) 
according to Eqs. ^ and ( [4c| . 

In Eqs. ([T]), the feedback term E}){t) has not been speci- 
fied yet. In the following, we introduce the feedback term 
such that it models a Fabry- Perot resonator. As opposed 
to the original LK model where only a single external mir- 
ror is considered, we take an external FP resonator with 
multiple reflections into account: 

Eb{t) (3) 

CO 

= Ke-'^ [^(^ -6-mr)-E{t-6-{m^ l)r)] 

= Ke-'"^ [E(t -S)-E{t-S- r)] + REb{t - r). 

K is the feedback strength, r is the delay time (cavity 
round trip time), is a memory parameter (mirror re- 
flectivity), 5 denotes the latency time originating from a 
single round trip between the laser and the resonator (see 
Fig. [T]), and cp is the feedback phase which results from 
the associated optical phase shift. 

The latency time ^, i.e., the propagation time between 
the laser and the FP, is correlated to the phase (p by the 
relation (f = Qq^^ where Qq is the frequency of the emit- 
ted light. However, we consider the two parameters cp 
and S as independent variables because the phase cp can 
be tuned by subwavelength changes of the separation be- 
tween laser and FP, on which scale the slowly varying am- 
plitude E^ which depends upon J, does not change. The 
effect of latency in time-delayed feedback was already 
studied in a general context in Refs. [22 , 24, 27, 28 . 

Throughout this work we use resonant feedback from 
the FP, i.e., the additional phase shift of the electric field 
by round trips in the FP resonator, denoted by ^, is 
assumed to be 27rm, where m is an integer. Otherwise the 
control term would be modified by an additional factor 
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e in the second term, and would not vanish for E{t 
S) = E{t-S-T). 



III. LINEARIZATION AROUND THE EASING 
FIXED POINT 



To investigate the stabihty of the cw laser emission, we 
will first consider the system without the feedback term. 
The uncontrolled system given by Eqs. ([T]) has a trivial 
fixed point at {n = I^E = 0). This fixed point is not of 
interest, because it describes a non-lasing state, i.e., the 
electric field is zero. 

Another fixed point, which describes a lasing state is 
located at {n = 0^ E = V^e^^). Now, the phase of the 
electric field can be arbitrarily fixed to = 0, such that 
_ (solitary laser mode). This leads to symmetry 
breaking of the rotational (SI) symmetry of Eqs. ([T]) with 
respect to complex E. 

Using the abbreviations 



E{t) = 

r = 



T [no + x{t) + iy{t)] ^ 



1 + / 1 



dk 
dn 



n=0 



(4a) 
(4b) 

(4c) 



with real- valued x and y, the fixed point is located at 
(n = 0,a; = 0,y = 0). A linearization around this point 
leads to the following system of equations: 



-4!^o 0\ ( n{t)\ 
x{t) . (5) 
y V y{t) J 



Note that an additional time scale transformation t 
(T /2)t was performed to eliminate the parameter T from 
the equations. Thus the rescaled time t is related to the 
physical time s by 





t 



s T 

Vr 2" 



s 
2^ 



(6) 



Also, the function k{n) is no longer present in the lin- 
earized equations. Instead, a parameter F was introduced 
that includes the differential gain dk/dn evaluated at the 
fixed point. The relation between T and the parameter 
no used in k{n) according to Eqs. ([2| and (4c) is shown 
in Fig. [2|b). The stability properties of the lasing fixed 
point are evaluated within this linearized model in the 
following. 

In order to investigate the stability of the fixed point 
in the uncontrolled system, we consider the characteristic 
equation for the eigenvalue A which is given by: 








Aid 



(7a) 



= A(A2 + 2FA + 41]^), (7b) 
where Id is the identity matrix. The solutions of this 
equation are given by 



A = 0, 



(8a) 
(8b) 



The Goldstone mode A = does not contain information 
about the stability of the fixed point and is only present 
due to the rotation symmetry. Therefore, only the second 
solution, given by Eq. (8b), is of interest. 



Note that the imaginary part y4^^o^^^' which cor- 
responds to the frequency of the relaxation oscillation, 
exists only for \2Qo\ > |F|. Under this condition the 
fixed point is a focus. Since — F is the real part of the 
eigenvalue A, the focus is stable for any F > and unsta- 
ble otherwise. Starting with an unstable focus (F < 0), 
we investigate the effects of the extended time-delayed 
feedback control on the stability of the fixed point in the 
following. 

Adding the control term Eq.([3| with properly rescaled 
K to the linearized equations leads to the following sys- 
tem: 





^m=0 





K \ cos (f sin (p 
- sin if cos cp 



x{t -S- mr) -x{t-S-{m^ l)r)] . (9) 



E^=o W -S-mr)-y{t-5-{m^ l)r)] 



4 



Using an exponential ansatz exp (At) for all three variables x, and n leads to the characteristic equation 



= (2r + A) 



Ke 



-A6 



1 - e 



-At 



+4ng (A^Ke 



-A5 



1 - Re-^^ 

1 - e-^^ 
1 - Re-^^ 



• + 2AKe 



-A6 



1-e 



-At 



COS cp + aKe 



-AS 



I -Re 

e 



1 .-A 



1-Re- 



- sm (p 



(10) 



In our simulations, we use the following parameters, 
which were chosen close to the values of Ref. [25^ : Qq = 
0.06, a = 5, r = -0.01 (corresponding to T = 500, 
/ = 1.8, A = 1, W = 0.02, no = -0.033933, and 
ko = 0.993075). Thus the intrinsic period of the un- 
controlled unstable focus is Tq ^ tt/Qq ^ 52. 

Note that the value of ko is determined by the other 
three parameters via the constraint k{0) = 1. With the 
values given here, the domain of control will be investi- 
gated in the {K^ r)- and (K, i?)-plane in the following. 

In the special case of = and ^ = the boundary 
of the control domain can be obtained analytically in 
a similar way as in Ref. [24 for the generic normal form 
model. The characteristic equation (10) can be factorized 
into 
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-At 



l-Re- 



(2r + A) A + 



1 



1 - Re-^^ 



(11a) 
-A^ll (lib) 



where Eq. (11a) corresponds to the Goldstone mo de in 
the uncontrolled system and therefore only Eq. (lib) 
is considered for the stability analysis. Separating 
Eq. (lib) into real and imaginary part and using the 
Hopf condition Re(A) = yields after some trigonometric 
manipulations the parametric representation of the con- 
trol domain boundaries in the (K^r) parameter plane, 
parametrized by the imaginary part Im(A) =: q of the 
complex eigenvalue: 



K{q) 
ri{q) 
^2(9) 



(1 + R) [q^ + lent + V (r^ - 2nl)] 



arcsin [a{q)] + 2m7r 



[a{q)] + (2m + 1) tt 



(12a) 
(12b) 
(12c) 



with 



a{q) = {K{q)q{R-l) 

[g2 + 4(r-i]o)(r + i^o)]}/ 
{(i^(^)^ + i?V) (g^ + 4r^) 

-SnlR {Rq^ - 2K{q)T) + WR^Q^} , 

where m is any nonnegative integer. This parameter had 
to be introduced due to the multiple leaves of the arcsin 
function. 



IV. NUMERICAL RESULTS 

The stability of the lasing fixed point is given by Eq. 
We solved this equation via Newton's method. 



(10). 



Since the transcendental equation has an infinite number 
of roots, we scanned the complex plane as initial condi- 
tions of the root-finding algorithm to locate the eigen- 
value with the largest real part. The parameter space 
consisting of r, R^ and S is five-dimensional for 
fixed r and Qq. To visualize the domain of control, 
we consider two-dimensional sections of this parameter 
space. 
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FIG. 3: (Color online) Domain of control according to 
Eq. ( |10[ ) in the {K, T)-plane for different values of cp. The 
greyscale (color code) denotes the largest real part Re (A) of 
the eigenvalues A, only negative values are plotted. Pan- 
els (a), (b), (c), and (d) correspond to cp = 0, 7r/16, 7r/8, 
and 37r/16, respectively. Other parameters: F = —0.01, 
Ho 0.06, a = 5, = 0.7, and ^ = 0. 

In Fig. [sj the domain of control is shown in the (i^, r)- 
plane for different values of the phase: (p = 0, 7r/16, 
7r/8, and Stt/W in panels (a), (b), (c), and (d), respec- 
tively. The greyscale (color code) denotes the largest 
real part of the eigenvalues and is therefore a measure of 
stability. Note that only values of Re (A) < are plot- 
ted, thus the shaded regions correspond to a stable lasing 
fixed point, i.e., a stable cw output. The control domains 
form tongues separated by (white) regions of no control 
around r = uTq with n integer, just as in the generic 
model studied in Refs. |^ El]. It can be seen that the 
domain of control shrinks with increasing phase. Here, 
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the domains of control are cut off by boundaries from 
the upper left and right for increasing phase, leading to 
overall smaller regions of stability. The tongues of stabi- 
lization are also slightly distorted towards smaller values 
of r. Additionally, in this picture, the regions of optimum 
stability, denoted by bright (yellow) color, are shifted to- 
wards larger values of the feedback gain K. 
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FIG. 5: (Color online) Domain of control in the (X, T)-plane 
for different values of 5 and fixed Lp — 0. The greyscale (color 
code) denotes the largest real part Re (A) of the eigenvalues 
A, only negative values are plotted. Panels (a), (b), (c), and 
(d) correspond to J = 0, 2.5, 5, and 7.5, respectively. Other 
parameters as in Fig. [3] 



FIG. 4: (Color online) Domain of control in the (i^^, r) -plane 
for different values of Lp. The greyscale (color code) denotes 
the largest real part Re (A) of the eigenvalues A, only negative 
values are plotted. Panels (a), (b), (c), and (d) correspond 
io if — 297r/16, 157r/8, 3l7r/16, and 27r, respectively. Other 
parameters as in Fig. [3] 

Figure [4] shows the domain of control in the {K^ r)- 
plane for values of (p in the range [I.Stt, 27r], which corre- 
sponds to negative phases [— 0.27r, 0] because the feedback 
of Eq. ([3| shows a 27r-periodicity in the phase variable. 
Note that the minimum feedback is only slightly changed 
by decrease of the phase. However, at larger values of K 
the control domain is cut off by an asymptotic boundary 
from the upper right, leading to a strong decrease of the 
maximum feedback gain with decreasing cp. In addition, 
the successive tongues of stabilization at larger values of r 
vanish altogether. Also, the region of optimum stability, 
denoted by bright (yellow) color shrinks for decreasing cp. 

Next, we will investigate the role of the latency time in 
the {K, r)-plane. In Fig. [sj the domain of control in the 
{K, r) -plane is depicted for different values of the latency 
time, i.e., ^ = 0, 2.5, 5, and 7.5, and fixed cp = 0. For 
larger latency times, the domains of control shrink and 
it can also be observed that they are bent down towards 
smaller values of the time delay r. Note that the regions 
of optimum stability, denoted by bright (yellow) color, 
are only slightly affected by the change of the latency 
time. 

All figures shown here were obtained for a fixed value 
of the memory parameter R = 0.7 as used in Ref. [4]. 
To investigate the dependence of the control on we 
display the domains of control in the {K, i?)-plane for 
different values of the phase {cp = 0, 7r/8, 7r/4, and Stt/S) 
and fixed time delay r = 26 in Fig. [6] This value of r 



was chosen based on the results from the generic model 
considered in Refs. 211 [24], where it was shown that the 
optimum time delay is given by r = To/2 = 7r/Im(Ao), 
where Aq denotes the eigenvalue of the uncontrolled sys- 
tem. In the LK model, the imaginary part of the eigen- 
values in the uncontrolled system is given by Eq. (8b), 
i.e., Im(Ao) = ^/4D^^^^^. This leads to an optimum 
time delay: 



'opt 



(13) 



which yields for our parameters Topt ^ 26. 

Now, in Fig. [6j it can be seen that the domain of con- 
trol in the (i^, i?) -plane has maximum size for = for 
this choice of the time delay r. [See panel Fig. [6]^ a).] 
For increasing phase, the domain of control shrinks while 
moving to the upper right. Stability is then only achieved 
in a small region at large values of K and R. [See panels 
Fig.gb) to (d).] 

For the generic model as used in Ref. [22] , it was shown 
that stability is enhanced in the case of nonzero phase if 
the delay time is chosen different from its optimum value 
of half the intrinsic period Tq. Therefore, we consider 
the domain of control in the (i^, i?)-plane for a value of 
r = 5 ^ O.lTo in Fig. [T] and for the same values of the 
phase cp as in Fig. [6] Here, the domain of control is al- 
ready smaller in the case cp = 0. Increasing the phase, the 
domain of control shrinks in i^-direction, but is enhanced 
slightly in the i?-direction, leading to an overall larger do- 
main of control than for the optimum r = 26. However, 
in the generic normal form model [24], the benefit of the 
non-optimum time delay was much more pronounced. 

Similar investigations as performed here for nonzero 
phase, are reported in the following for nonzero latency 
time S and fixed (p = 0. In Fig. |8| the domain of control 
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FIG. 6: (Color online) Domain of control in the (K, i?)-plane 
for different values of cp and fixed optimum time delay r = 26. 
The greyscale (color code) denotes the largest real part Re (A) 
of the eigenvalues A, only negative values are plotted. Panels 
(a), (b), (c), and (d) correspond to cp = 0, 7r/16, tt/S, and 
37r/16, respectively. Other parameters as in Fig. [s] 
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FIG. 7: (Color online) Domain of control in the {K, i?)-plane 
for different values of cp and fixed time delay r = 5. The 
greyscale (color code) denotes the largest real part Re (A) of 
the eigenvalues A, only negative values are plotted. Panels 
(a), (b), (c), and (d) correspond to cp = 0, 7r/16, 7r/8, and 
Stt/IG, respectively. Other parameters as in Fig. [s] 



in the (i^, i?)-plane is shown for different values of the la- 
tency time and fixed r = 26. For better comparison, we 
choose the same latency times as in Fig. [5] It can be ob- 
served that the domain of control shrinks for increasing 5. 
The lower right boundary is shifted upwards. Therefore, 
control is only possible for large values of if 5 is large. 
The minimum feedback gain is also increased, since the 
left boundary is shifted to the right. The regions of op- 
timum stability, denoted by bright (yellow) color shrink 
for increasing 5^ leading to a deterioration of the control. 

Previously, it was shown in the generic model of an 
unstable focus |24j that proper tuning of the latency can 



FIG. 8: (Color online) Domain of control in the (K, i?)-plane 
for different values of 6. The time delay and feedback phase 
are fixed at optimum r = 26 and cp = 0, respectively. The 
greyscale (color code) denotes the largest real part Re(A) of 
the eigenvalues A, only negative values are plotted. Panels 
(a), (b), (c), and (d) correspond to ^ = 0, 2.5, 5, and 7.5, 
respectively. Other parameters as in Fig. [3] 
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FIG. 9: (Color online) Domain of control in the (K, i?)-plane 
for different values of 6. The time delay and feedback phase 
are fixed at r = 5 and cp = 0, respectively. The greyscale 
(color code) denotes the largest real part Re(A) of the eigen- 
values A, only negative values are plotted. Panels (a), (b), 
(c), and (d) correspond to ^ = 0, 2.5, 5, and 7.5, respectively. 
Other parameters as in Fig. [3] 



compensate for a bad choice of the time delay r. There- 
fore, we depict the domain of control in the (i^, i?)-plane 
in Fig.[9]for a value of r = 5 and the same values of 5 as in 
Fig. [Sj Note that the domain of control is smaller than in 
Fig. [8] for ^ = 0. However, the size of the domain of con- 
trol is greatly enhanced for increasing latency time. The 
region of stability is bent towards smaller values of the 
memory parameter allowing for control at a smaller 
value of the feedback gain for given R. The region of op- 
timum stability is located at a large value of R for S = 0. 
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This region also moves towards lower values of R and 
grows for increasing 5. 




FIG. 10: (Color online) Domain of control in the (K, (p)-p\ane 
for different values of r and fixed 6 = 0. The grey scale (color 
code) denotes the largest real part Re (A) of the eigenvalues 
A, only negative values are plotted. Panels (a) and (b) corre- 
spond to r = 26 and 5, respectively. Other parameters as in 
Fig.H 

To investigate the dependence of the domain of con- 
trol on the choice of the phase cp further, we consider an- 
other two-dimensional projection of the five-dimensional 
control-parameter space parameterized by feedback gain 
K and the feedback phase ip. This section is depicted 



in Fig. 10 for two different values of the time delay and 
fixed 6 = 0. In panel (a), the time delay is chosen as 26, 
which is the optimum r according to Eq. (13). Here, it 
can be seen that the optimum phase is located at slightly 
negative values for small values of the feedback gain up to 
K ^ 0.05. Increasing the optimum phase changes its 
sign and is now located at small positive values of (p. For 
the case of r = 5, which is depicted in panel (b), stability 
is overall enhanced drastically. The bright (yellow) areas, 
denoting regions of optimum stability, are located at neg- 
ative (f for small K up to K ^ 0.1. Control is possible 
even for a small value of (f below — 7r/2, if the feedback 
gain is tuned exactly to the small range of K ^ 0.01. For 
larger feedback gain with K > 0.1, the optimum value of 
(f is located at positive values. The region of optimum 
stability is located at large values of K around 0.2. The 



shape of the control domain in Fig. 10 is markedly dif- 
ferent from that in the generic normal form model (see 
Figs. 6(c) and 7(c) in Ref. [2T), but appears to be in 
line with full device simulations within a travelling wave 
model [5]. 



V. BANDPASS FILTERING OF THE 
FEEDBACK SIGNAL 

A bandpass filter is widely used in experiments with 
optical systems, as well as in theoretical treatments of the 
LK model [Ml [29] . It is used to suppress unwanted fre- 
quencies in the feedback loop but can also fundamentally 
influence the dynamics of the laser system [30 . There- 
fore it is useful to include a corresponding term into 
the theory. Experimentally, the filter is realized by a 
Fabry-Perot interferometer, which can approximatively 
be modelled by a Lorentzian. Therefore, the bandpass 



filter is introduced by using the transfer function T{lj) 
of a Lorentzian in Fourier space. The transfer function 
acting on the electric field E can be written as 



1 



7 



(14) 



where coq denotes the peak of the transfer function and 
7 the full width at half maximum. In Fourier space, the 
bandpass filter alters the complex electric field E{t) as 
follows 



T{u)E{ij) 
1 

7 



E{lo) 



(15) 



with the filtered electric field E. Transforming back from 
Fourier space yields a differential equation for E(t): 



dE{t) 
dt 



= (ic^o-7)^W+7^W, 



which is added to the model equations ([T]), 

dE T , ^ ^ ^ . X 

(1 + ia) nE-Et{t), 

- [l^n]k{n)\Ef 
= (ic^o-7)^W+7^W, 



dt 
dn 

'dt 
dE 
~dt 



2 



(16) 

(17a) 
(17b) 
(17c) 



while the filtered field E{t) is used instead of E{t) in the 
feedback term ([3|: 



Eb{t) 



Ke-'"^ J2 [E{t -S-mr)- E{t 



{IS) 

^-(m + l)r)] 



Similar to the unfiltered system, the complex filtered 
variable E(t) can be split into real and imaginary parts 
denoted by x{t) and y{t), respectively. This leads to addi- 
tional differential equations for the filtered variables x{t) 
and y{t): 



dx(t) 
dt 

dm 

dt 



-f{x{t)-x{t))-^oy{t). (19a) 
j{y{t)-y{t))^^ox{t). (19b) 



where we rescaled the parameters 7 and loq by 2/T. The 
time has also been rescaled as in Eq. ([5|. 

For the sake of simplicity, we restrict ourselves to the 
case of cp = and (5 = in the linear stability analysis, 
which corresponds to the case of zero phase and no ad- 



ditional latency. Using Eqs. [T9| for the filtered variables 
X and y and the linearization from Eq. ([s]), the charac- 
teristic equation of the controlled system with bandpass 
filter can be obtained: 
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-Ar 



= 



Aa-fuJonlK- -—- + KujI {2TK + + ml 

1 — tie 

1 - ' 



-Ar 



I -Re 



:^+A(7 + A) 



■ (2r + A) 



l-Re- 



^+A(7 + A) 



(7 + A) 



(20) 



The eigenvalues A, which are the roots of the charac- 
teristic equation, are calculated numerically. Details on 
the numerical procedure are given in the beginning of 



frequency components are cut off and thus, more infor- 
mation of the system is lost due to the filter. 



Sec. IV The parameter space is five-dimensional and pa- 
rameterized by r, i?, 7, and uoq. Results are shown 
in projections onto the {K^r) parameter plane in the fol- 
lowing. 
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FIG. 11: (Color online) Domain of control in the (K, T)-plane 
according to Eq. ( [2Q| ) for different values of 7 = 100, 1, 0.1, 
and 0.05 in panels (a), (b), (c), and (d), respectively. a;o = is 
fixed. The greyscale (color code) denotes the largest real part 
Re (A) of the eigenvalues A, only negative values are plotted. 
Other parameters as in Fig. [S] 



In Fig. 11 , the domain of control is shown in the (i^, r)- 
plane for fixed coq = and different values of the filter 
width: 7 = 100, 1, 0.1, and 0.05. The choice of coq corre- 
sponds to the case of the low-pass filter because the peak 
of the transfer function is shifted to zero. The greyscale 
(color code) corresponds to the largest real part of the 
eigenvalues. Only negative values are displayed. For 
large values of 7, the picture is almost identical to the 
case of the unfiltered system shown in panel (a) of Fig. [s] 
This can easily be understood because, for large cut-off 
frequency 7, the electric field passes the filter almost un- 
changed. Decreasing the filter width to 7 = 1 has only 
little effect on the shape of the domain of control, but the 
region of best stability, denoted by bright (yellow) color, 
becomes larger. Decreasing 7 further, the domain of con- 
trol shrinks and is bent down slightly towards smaller 
values of r. This suggests that the optimum r changes, 
which will be investigated later in the (i^, i?)-plane. Fur- 
ther increase of 7 leads to smaller domains because higher 
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FIG. 12: (Color online) Domain of control in the (K, T)-plane 
according to Eq. ( [20| for different values of ooo with fixed 7 = 
0.1. The greyscale (color code) denotes the largest real part 
Re (A) of the eigenvalues A, only negative values are plotted. 
Parameters: F = -0.01, Ho = 0.06, a = 5, and i? = 0.7. 
Panels (a), (b), (c), and (d) correspond to c^o = 0, 0.001, 
0.005, and 0.01, respectively. 



Figure 12 depicts the domain of control in the {K^ r)- 
plane for different values of the filter center frequency 
ujo = 0, 0.001, 0.005, and 0.01 in panels (a), (b), (c), and 
(d), respectively. The value of 7 is fixed at 7 = 0.1. The 
increase of coq has only little effect on the first tongue of 
stability. However, the second tongue is cut off from the 
upper right with larger ujq^ leading to a smaller range of 
possible values for the feedback gain K. This tongue gets 
also slightly thicker in r direction for increasing coq. 



VI. DELAY-INDUCED MULTISTABILITY 

In all preceding considerations, the effect of time- 
delayed feedback control on the stability of the lasing 
fixed point was investigated only locally, since the sys- 
tem was linearized around the fixed point. Results from 
dynamical simulations of the full nonlinear system, given 
by Eqs. ([T]), with feedback according to Eq. (|3|, are re- 
ported in this Section. 

Besides the fixed point associated with cw emission, 
the model without control, described by Eqs. ([T]), has a 
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LU 




FIG. 13: Absolute value of the electric field vs reduced carrier 
density: Limit cycle of the undamped relaxation oscillations 
in the modified LK model according to Eq. ([T]) and Eq. ^ 
without feedback. Parameters: T = 500, a = 5, / = 1.8, 
A = 1,W = 0.02, and no = -0.034. The arrow indicates the 
direction of the trajectory. The cross denotes the lasing (cw) 
fixed point. 
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limit cycle in the {\E\^ n)-plane, which corresponds to un- 
damped relaxation oscillations (or intensity pulsations). 
The amplitude of this limit cycle, which is born by a 
Hopf bifurcation at F = 0, depends on the choice of the 
parameters. For the same parameters as in Sec. |III[ the 
limit cycle is depicted in Fig. [Tsj The unstable fixed 
point is located at (n = 0, |£^| = 1.342) for this choice of 
parameters. 

In an experiment, the case may be of interest where 
the feedback is activated when the laser is operating 
at this limit cycle of relaxation oscillations. Therefore, 
for numerical simulations we choose initial conditions on 
this limit cycle: E{t = 0) = -0.731767 - i 0.016891, 
which corresponds to x{0) = —0.0927256 and y{0) = 
—0.00075538. The initial values of the reduced carrier 
density is chosen as n(0) = —0.037433. 

In the presence of a feedback term, the system equa- 
tions become a set of delay differential equations. There- 
fore, the initial conditions must be specified for the time 
interval [— r, 0). For this interval, we choose the values of 
the limit cycle shown in Fig. [13] Depending on the feed- 
back gain K and the time delay r, the system exhibits di- 
verse scenarios under the influence of feedback according 
to Eq. (|3|. Some trajectories of the system in the n)- 
plane are displayed in Fig. 14 for fixed time delay r = 26 
and different values of the feedback gain K. The memory 
parameter is fixed to = 0.7. A variety of dynamic sce- 
narios ranging from limit cycles to chaotic attractors can 
be observed depending on the choice of K. Panel (a), in 
the absence of control {K = 0), shows the limit cycle of 
undamped relaxation oscillations as Fig. [13] Increasing 
the feedback gain to K = 0.03 (b), a limit cycle of period 
2 evolves. For a value of = 0.06 (c), chaotic behavior is 



FIG. 14: Dynamics with feedback for different values of the 
feedback gain K. Initial conditions are chosen on the limit 
cycle of relaxation oscillations (Fig. |13|). The time delay and 
memory parameter are fixed at r = 26 and R = 0.7, respec- 
tively. = 0, ^ = 0. Panels (a), (b), (c), (d), (e), and (f) 
correspond to i^^ = 0, 0.03, 0.06, 0.1, 0.12, and 0.15, respec- 
tively. The rotational direction of the trajectories is the same 
as in (a) for all other panels. Other parameters as in Fig. 13 
Note that the trajectories are shown only for a limited span 
of time. 



observed. At = 0.1 (d) a small limit cycle appears, but 
this limit cycle does not oscillate around the uncontrolled 
fixed point, which is located at n = 0. Upon further in- 
crease the limit cycle shrinks (e), and finally disappears 
in an inverse Hopf bifurcation from a different fixed point. 
A different limit cycle is observed for K = 0.15 (f). For 
other parameters, a variety of complex scenarios is found. 
This is consistent with recent findings by Tronciu et at. 
^, who demonstrated rich multistable dynamic scenar- 
ios of feedback-induced stationary external cavity modes 
(rotating waves E{t) = Ege^^^'^ and further Hopf bifurca- 
tions of those) if a certain critical feedback strength Kc is 
exceeded. Note that the lasing fixed point (solitary laser 
mode) at n = is stable for all values of > in Fig. 14 
thus there exists multistability between cw laser emis- 
sion and (periodic or chaotic) intensity pulsations, and 
it depends upon the initial conditions which attractor is 
asymptotically reached. If the initial condition is chosen 
in the vicinity of the fixed point, where the linearization 
is valid, time-delayed feedback control still works to sta- 
bilize cw emission. This can be achieved, for instance, by 
first operating the laser below the Hopf bifurcation and 
then gradually increasing the pump current, and hence 
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r, as will be shown in the next Section. Note that delay- 
induced multistability is common in many systems [31^ . 



VII. TURN-ON DYNAMICS 

The results of the preceding Section might suggest that 
stabilization of the lasing fixed point (cw emission) in the 
full nonlinear system is difficult if not impossible for the 
majority of choices for the control parameters r, 
5^ 0, 7, and uoq. Driving the system into the lasing fixed 
point is only possible for initial conditions in a vicinity 
of the lasing fixed point. An experimental realization 
of such initial conditions will most probably be compli- 
cated. Choosing other initial conditions, for example on 
the limit cycle of relaxation oscillations, as shown in the 
previous Section, leads to unwanted multistability in the 
presence of control. 

We have therefore simulated a turn-on scenario which 
is experimentally both realistic and feasible. As a result 
of the slowly increasing instability the feedback is able to 
keep the dynamics in the very vicinity of the lasing fixed 
point. Even with a small amount of noise this behavior 
is expected to persist, since the basin of attraction of the 
fixed point is finite. 
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FIG. 15: (Color online) Dynamics of the system during a turn- 
on process without feedback (iC = 0). The dashed (red), solid 
(green), and dotted (blue) curves correspond to the injection 
current /, the carrier inversion n, and the absolute value of 
the electric field \E\^ respectively, vs rescaled dimensionless 
time. Other parameters as in Fig. [3] 

The turn-on scenario is depicted in Fig. [15] for the 
case of no control. The laser is first operated below 
the laser threshold with the injection current chosen as 
/ = —1. With initial conditions chosen arbitrarily as 
X = y = 10~^ and n = the laser relaxes into the non- 
lasing fixed point. The injection current is then linearly 
increased until the value of / = 1.8 is reached, which was 
considered in the previous Sections of this work. Pass- 
ing the laser threshold at / = 0, non-lasing and lasing 
fixed points change stability, thus the system swaps to 



the lasing fixed point with n = and nonzero E. In- 
creasing the current further beyond the value of / = 0.3, 
the lasing fixed point becomes again unstable due to the 
internal dispersive feedback (see Eq. (|8b[)). Reaching 
/ = 1.8, the system resides in the limit cycle shown in 
Fig.[l3| The full turn-on process from / = —1 to 1.8 takes 
9 X 10 in the rescaled time units according to Eq. (|6| 
and used in these simulations. Assuming a photon life- 
time of the order of ^ 10~^^s, this would relate to a 
physical time of 5turn-on ^ 1.8ms for the turn-on ramp 
using s = {2Tc/T)t = 2rpt, where t is the rescaled dimen- 
sionless time and s is the time with physical units. Note 
that Tc = 5 X 10~^s for the chosen values of T and r^. 
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FIG. 16: (Color online) Dynamics of the system during a 
turn-on process with feedback {K = 0.1, r = 26). The 
dashed (red), solid (green), and dotted (blue) curves corre- 
spond to the injection current /, the carrier inversion n, and 
the absolute value of the electric field |^|, respectively. Other 
parameters: R — 0.7, = 0, and ^ = 0. 

We now consider the system in the presence of control. 
For experimental feasibility the feedback is active during 
the full turn-on process. We consider an exemplary case 
with r = 26, K = 0.1, and R = 0.7. Phase and latency 
are chosen as zero, the bandpass filter is turned off. The 
dynamics during the turn-on process is shown in Fig. 16 
It can be seen that the system remains in the vicinity 
of the lasing fixed point even beyond the point / = 0.3. 
This is achieved since the feedback is already active as 
the system passes this onset of instability. As the system 
is thereby always kept near the fixed point the results of 



the linear stability analysis of Section III can be applied 
here. Another evidence for the success of the control is 
the fact that after the turn-on process the control force 
vanishes. In the simulations of Section |VI[ where only 
delay-induced orbits were created, the control force did 
not vanish. 



In Fig.[T7| the analytic results from Section III are plot- 
ted together with the results from turn-on simulations as 
shown before. In the lower part of the picture for rela- 
tively small values of r, the results from the simulations 
match the results from the linear stability analysis. How- 
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FIG. 17 : (Color online) Domain of control according to 
Eq. |Tq| in the {K, T)-plane compared to results from turn-on 
simulations. The greyscale (color code) denotes the largest 
real part Re (A) of the eigenvalues A, only negative values are 
plotted. Circles show results from simulations. Full (blue) cir- 
cles correspond to successful control of the fixed point, empty 
(red) circles denote no control. Parameters: a = 5, = 0.7, 
= 0, and 6 = 0. 

ever, in the upper part of the lower tongue of control and 
in the whole upper tongue of control, the simulations 
are not successful. The large value of the time delay 
makes control less effective at the onset of instability of 
the fixed point (at / = 0.3). This drives the trajectory 
slightly away from the fixed point, which can afterwards 
not be corrected by the control. Choosing a slower ramp 
for the turn-on process retains control for larger values 
of r. We have performed exemplary simulations with a 
turn-on time ten times longer than before (now 9 x 10^ 
in rescaled dimensionless units) using K = 0.03 for a dif- 
ferent value of the time delay. For r = 35, control was 



then successful. 



VIII. CONCLUSION 

In conclusion, we have shown that time-delayed feed- 
back provides a valuable tool for the suppression of un- 
wanted intensity pulsations in a semiconductor laser, 
which can conveniently be realized by a Fabry-Perot res- 
onator. We stress that the efficiency of this method has 
already been demonstrated experimentally [4]. We have 
discussed the effects of the control scheme in the frame- 
work of a modified Lang-Kobayashi model. Our results 
were obtained by a linear stability analysis of the fixed 
point of the delayed system. As a modification of the 
original controller design, we have taken into account 
an additional control loop latency as well as a variable 
phase-dependent coupling. We have shown how these 
extensions affect the domain of control in various pro- 
jections of the parameter space. Futhermore, we have 
investigated the effects of a bandpass filter added in the 
feedback. In addition, we have presented simulations of 
the full nonlinear system, which may exhibit multistabil- 
ity, so that simple feedback control fails. In this case, 
control may be achieved only by slowly increasing the 
injection current until the desired value is reached. This 
corresponds experimentally to turn-on of the laser with 
a slow injection current ramp. 
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